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Abstract 

In the computation of a Grobner basis using Buchberger's algorithm, a key issue for 
improving the efficiency is to produce techniques for avoiding as many unnecessary 
critical pairs as possible. A good solution would be to avoid all non-minimal critical 
pairs, and hence to process only a minimal set of generators of the module generated 
by the critical syzygies. In this paper we show how to obtain that desired solution 
in the homogeneous case while retaining the same efficiency as with the classical 
implementation. As a consequence, we get a new Optimized Buchberger Algorithm. 
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1 Introduction 



Ever since practical implementations of Buchberger's famo us algorithm for 



computing Grobner bases became feasible (jBuchbergerl . ll965f b it has been clear 
that, in order to improve the efficiency of this a lgorithm, on e needs to avoid the 
treatment of as many critical pairs as possible. Ilhiclilx'v.'cil (ll 97flh studied this 



probl em for the first time, and later in (|Buchbergerl . 119851 ) and ((Gebauer and Moller 



1987 ) his results were sub stantially improved and expanded. Nevertheless, 



Gebauer and Mollerl (|1987f ) showed that their method did not always produce 
a minimal set of generators of the module generated by the critical syzygies. 
However, their method was very efficient and yielded an almost minimal set 
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of critical pairs. Since then, many kinds of optimizations of Buchberger's al- 
gorithm have been found, in particular by implementers of computer algebra 
systems. But the problem of efficiently minimalizing the critical pairs has gone 
largely unnoticed and seems to be overdue for a solution. Indeed, that is the 
main objective of this paper. 

To achieve our goal, we proceed as follows. Foremost, we need a detailed un- 
derstanding of the entire process of computing Grobner bases, in particular in 
the homogeneous case. An algorithm for simultaneously computing a Grobner 
basis and a minimal system of generators contained in it is fine-tuned when 
the input is a reduced Grobner basis. Then this result is applied to critical 
syz ygies, using the fact that w e show how the old criteria M(i,j) and F(i,j) 



of ([Gebauer and Mollerl . 119871 ) yield a reduced Grobner basis of the module of 
syzygies of the leading terms. Besides, when applied to this special case, the 
algorithm admits many subtle optimizations. In the end, we really achieve the 
goal of minimalizing the critical pairs efficiently. 

Now, why do we think that what we achieved is important? 

The first reason is theoretical curiosity. It is common knowledge among the im- 
plementers of Buchberger's algorithm that the criteria of Gebauer and Moller 
almost produce a minimal set of critical pairs. We wanted to see whether that 
vox populi is really true. Of course one could use a standard minimalization 
process to produce minimal sets of critical pairs, but this method could only 
handle small examples. Instead, we observed that, after applying two of the 
criteria of Gebauer and Moller, a reduced Grobner basis of the module of syzy- 
gies of the leading terms is obtained. Then we were able to see the difference 
between the reduced Grobner basis and a minimal set of generators of this 
module, and how this difference depends on the size of the example. 

Another important reason is that we wanted to be able to compute a minimal 
set of generators of this module with the same efficiency as the usual applica- 
tion of the Gebauer-Moller criteria. And we wanted to do it while computing 
a Grobner basis, so that we can replace the Gebauer-Moller criteria by our 
procedure. As we show in the last sections, we achieved this goal. 

A third reason is that our results hold in full generality, namely for Grobner 
bases of modules over positively (multi-) graded rings. Other optimizations 
of Buchberger' s algorithm, e.g. ideas using trivial syzygies (see for instance 
Faueer 3 £3)), do not hold in this generality. Moreover, we would like to 



point out that the pairs we discard are truly useless, whereas pairs between 
elements in a reduced Grobner bases which reduce to zero can still be useful 
for the computation of syzygies. 

Finally, the readers should know that the basic terminology is taken from the 
book of the second and third authors ( Kreuzer and Robbianol 2Q00h . 
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2 Some Background Material 



Since we are interested in optimizing Buchberger's algorithm in the homo- 
geneous case, we start by saying which gradings we consider. From now on 
let K be a field and P = K[x±, . . . , x n ] a polynomial ring over K . Moreover, 
let W G Mat TO>n (Z) be an mxn-matrix with integer entries. Then there exists 
exactly one Z m -grading on P such that every term t = x" 1 ■ ■ ■ x" n is homo- 
geneous of degree deg w (t) = W ■ (cti, . . . , a n ) tr . We say that P is (multi-) 
graded by W . The matrix W is called the degree matrix and its rows are 
called the weight vectors. 

For instance, the grading on P given by W = (1, . . . , 1) is the standard 
grading. For every d G Z m , the homogeneous component of degree d of P 
is Pw,d = ®dcg w (t)=d K ■ t. Given 5i, . . . ,5 r G Z m , the graded free P-module 
F = (B r i= i P(—8i) inherits a Z m -grading from P in the natural way. Again we 
say that F is graded by W . 

In order to be able to use these gradings in our algorithms, we need some 
positivity assumptions. 



Definition 1 Let P be graded by W, and let wi, . . . , w m be the rows of W . 

a) The grading given by W is called weakly positive if there exist integers 
ai, . . . , a m such that a\W\ + • ■ ■ + a m w m has all entries strictly positive. 

b) The grading given by W is called positive if rk(W) = m, if no column 
of W is zero, and if the first non-zero entry in each column of W is 
positive. 



Proposition 2 Let P be weakly positively graded by W, and let M be a 

finitely generated graded P-module. 

a) We have Pw,o = K and dim^ {My/4) < 00 for every d G Z' m . 

b) The graded version of Nakayama's lemma holds: homogeneous elements 
Vi, . . . , v s G M generate the module M if and only if their residue classes 
vi, . . . , v s generate the K -vector space Mj (xi, . . . , x n )M . In particular, 
every homogeneous system of generators of M contains a minimal one, 
and all irredundant homogeneous systems of generators of M have the 
same number of elements which is denoted by fi(M) . 



The proof of this proposition uses standard compute r algebra methods and is 
contained in ( Kreuzer and Robbiand . in preparation ). For practical computa- 
tions we need the somewhat stronger notion of a positive grading. The use- 
fulness of positive gradings is illustrated by the following characterizations. 
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Recall that a module ordering a on the set of terms T n (ei, . . . , e r ) of the 
graded free module F is called degree compatible or compatible with 

deg^y if the inequality deg w (tei) > Lex deg w (t'ej) implies te.{ > a t'tj for all 
t,t' G T n and all i,j G {l,...,r}. 

Proposition 3 Let P be graded by W , where W has Z -linearly independent 
rows and non-zero columns. Then the following conditions are equivalent. 

a) The grading on P given by W is positive. 

b) The restriction of Lex to the monoid Y = {d G Z m | Pw,d 7^ 0} is a 
well- ordering, i.e. every non-empty subset of V has a minimal element 
with respect to Lex. 

c) The restriction of Lex to the monoid T = {d G Z m | Pw,d 7^ 0} is a term 
ordering, i. e. every element d G V satisfies d > Lex . 

d) There exists a term ordering r on T™ which is compatible with deg w . 

e) There exists a module term ordering o on T n (ei, .. . ,e r ) which is com- 
patible with the grading given by W . 



Again we refer to ( Kreuzer and Robbianol in preparation ) for a proof of this 



proposition. As a consequence, it follows that positive gradings are weakly 
positive. Moreover, in a positively graded setting, we can prove the fmiteness 
of various algorithms in the usual way, i.e. by using the fact that there is no 
infinite sequence of homogeneous elements of strictly decreasing degrees. 

In the remaining part of this section, we use truncated Grobner bases to prove 
two very important technical tools, namely Corollary 8 and Corollary 10. We 
shall from now on assume that P is positively graded by W G Mat mi „(Z). 
Moreover, we let 8x,...,8 r G Z m , we let M be a finitely generated graded 
sub module of the graded free P-module F = ©[ =1 P(—Si) , and we let a be a 
module term ordering on T n (ei, . . . , e r ) , the set of terms in F . 

The following notation will turn out to be convenient. Given a subset S of 
a graded P-module and d G Z m , we let S<d = {v G S \ v homogeneous, 
do 

Sw( v ) ^Lex d} and S d = {v e S \ v homogeneous, deg w (v) = d} . 



Definition 4 Assume that G = {gi, . . . , g s } is a homogeneous cx-Grobner 
basis of M, and let d G Z m . Then the set G<d is called a c/-truncated 
Grobner basis of M, or a Grobner basis of M which has been truncated 
in degree d. 

For truncated Grobner bases, we now prove a characterization which is anal- 
ogous to the Buchberger criterion in the usual case. To this end, we need to 
explain what we mean by critical pairs and critical syzygies. 
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Given homogeneous elements gi, . . . , g s G M\ {0}, we let di = deg w (gi) for 
i = 1, . . . , s, and we let F' be the graded free P-module ©f =1 P(—di) . The 
canonical basis of F' will be denoted by {e 1; . . . ,e s }. Notice that we have 
deg w (Ei) = di for i = 1, . . . , s. Moreover, we write LM a (gi) = crfie^, where 
Cj G K \ {0} , where ^ 6 P, and where 7» G {1, . . . , r} . 



Definition 5 A pair G {l,...,s} such that 1 < % < j < s and 

7i = 7j is called a critical pair of (gi, . . . , g s ). The set of all critical pairs 
of (<?,.. • , <7s) is denoted by B. For every critical pair G B, the element 
a ii = 1Cm e A tj) £ » ~ Wy.?) gj . i s a syzygy of the pair (LM ff fe),LM ff (ft)). It is 
called the critical syzygy associated to the critical pair (i,j) . The set of all 
critical syzygies is denoted by E. 



Clearly, a critical syzygy <Jij is a homogeneous element of F' whose degree is 
precisely deg w (<Jij) = deg w (\cm(ti,tj)) + 5 yi . This degree equals the degree 
of the corresponding S- vector Sa = lcm ( f " f j) „. _ lcm ^ t » g. i n p. 

For every critical pair G B, we call deg w (<Jij) the degree of the critical 
pair. Then it makes sense to consider the set B<^ for every given d G Z m , and 
we observe that deg w ((Ty) > Lex ma.x{di,dj} for all G B. Finally, we re- 
mind the reader that NR^ c^V) denotes normal remainde r, i.e. the result of the 
division algorithm, as in ( Kreuzer and Robbianol . l200(l ). Definition 1.6.7. At 



this point, we are ready to formulate and prove the following characterization 
of truncated Grobner bases. 

Proposition 6 (Characterization of Truncated Grobner Bases) 

Let P be positively graded by W G Mat mjn (Z), let G = {gi, . . . , g s } be a 
set of non-zero homogeneous vectors which generates a graded submodule M 
°f ©[=i P(— 5i) ; an d ^ d G Z m . Then the following conditions are equivalent. 

a) The set G<d is a d -truncated a -Grobner basis of M. 

b) For every homogeneous element v G M<^ \ {0}, we have the relation 
LT,(v) E (LTM \ g e G< d ) . 

c) For all pairs G B<^ ; we have NR a g <d (Sij) = 0, where Q<d is the 
tuple obtained from Q = (gi,...,g s ) by deleting the elements of degree 
greater than d. 

Proof. Without loss of generality, we may assume that G<d = {gi, . . . , g s >} for 
some s' < s. It is clear that a) implies both b) and c). Now we show that b) 
implies a). By the assumption, we can find terms t' s , +1 , . . . , t' s „ of degree greater 
than d such that the set {LT CT (gi), . . . , LT a (g s >)} U {t' s / +1 , . . . , t' s „} is a system 
of generators of LT CT (M). We choose homogeneous elements /v+i, . . . , h s » 
in M such that LT^ (hi) = for i = s' + 1, . . . , s" . Then the set {gi, . . . , g s >, 
h s ' + i, . . . , h s /'} is a homogeneous cr-Grobner basis of M with truncation G<d- 
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It remains to prove that c) implies b). Let v G M<^ be a homogeneous 
non-zero element. Since {g x , ■ ■ ■ ,g s >} generates (M<d) , we can represent v as 
v = E|li fi9i, where ft is homogeneous of degree deg w (v) -deg w (gi) < Lex d. 
In order to prove LT CT (f ) G (LT^ffx), . . . , LT ff (q,,/)) , it is enough to proceed as 
in the proof of Proposition 2.3.12 of (|Kreuzer and RobbianoL 2000[ ). replacing 



Q by Q< d . □ 
This characterization has several useful applications. 

Corollary 7 Let G = {gi, . . . ,g s } be a homogeneous a -Grobner basis of the 
module M , and let d G Z m . Then G<d is a d -truncated a -Grobner basis of 
the module (M< d ) . 

Proof. Since G is a set of generators of M , the set G<d generates the module 
(M<d) . From Buchberger's Criterion we know that NR . ) g(5'y) = 0, for all 
pairs G B. If we have deg w (Sij) <Lex d here, the elements of G involved 
in the reduction steps SV,- — > all have degrees less than or equal to d. Hence 
we see that NR at g <d (Sij) = 0, and the proposition yields the claim. □ 

Corollary 8 Let d G Z m , let the elements of the tuple Q = (gx,...,g s ) 
form a d -truncated a -Grobner basis of M, and let g s+ \ G F be a homo- 
geneous element of degree d such that UY a (g s+ i) £ (LT a (gi), . . . , LT a (g s )) . 
Then {g±, . . . , g s+ i} is a d -truncated Grobner basis of M + (g s +i) ■ 

Proof. In order to prove the claim, we check condition c) of the proposition. 
For 1 < i < j < s such that deg w (Sij) < Lex d, we have NK a> g(Sij) = 
by the assumption and by Proposition 6. For i G {l,...,s} such that 
deg w (Si s+ i) = d, the fact that the pair (i,s + 1) has degree d implies that 
LT a (g s+ i) is a multiple of LT a (gi), in contradiction to the hypothesis. □ 

In the last part of this section, we prove an analogue of the preceding corol- 
lary for minimal generators. Recall that Proposition 2.b guarantees that all 
minimal systems of generators have the same length in the positively graded 
situation. 

Proposition 9 Let P be positively graded by W G Mat mjn (Z) ; let M be a 
graded P -module generated by homogeneous elements {g±, .. .,g s } , and as- 
sume that deg w {g x ) < Lex deg w (g 2 ) < Le x ■ • ■ <Lex deg w (g s ) . 

a) The set {gi, . . . , g s } is a minimal system of generators of M if and only 
if we have gi {g x , • • • , ft-i) for i = 1, . . . , s . 

b) The set {gi | i G {1, . . . , s}, gi (g x , . . . , ft-i)} is a minimal system of 
generators of M . 

Proof. First we prove a). If {g x , . . . , g s } is a minimal set of generators of M , 
then no relation of type g^ G (g x , . . . , ft-i) holds, since otherwise we would 
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have M = {g h g^, g i+1 , ...,g s ). Conversely, if {g u ...,g s } is not a mini- 
mal set of generators of M , then there exists an index % g {1 s} such that 

Qi g (ffi, . . . , ft-i, ft+i, • • • , g s ) ■ Using Corollary 1.7.11 of (|fCreuzer and Robbianol . 
2000), we obtain a representation g^ = J2j^ifj9j, where fj g P is homoge- 
neous of degree deg w (gi) - deg w (gj) for j g {1, . . . , % - 1, i + 1, . . . , s}. 

Since deg w (/j) > Lex for fj ^ 0, we see that deg w (gi) < Lex deg w (gj) im- 
plies fj = 0. Thus there are two possibilities. Either we have deg w (gi) > Lex 
deg w (gj) for all j such that fj ^ or there exist some indices j such that 
deg w (gj) = deg w (gi). In the first case, those indices j satisfy j < i by the 
assumption that the multidegrees of gi, . . . , g s are ordered increasingly, and 
therefore we get jj e (ft,..., 9i-i) ■ in the second case, the fj corresponding to 
those indices j are in K\ {0} . Let j max = max{j g {1, . . . , s} \ fj g if \ {0}} . 
We get the relation gj max g (pi, . . . , <7j max -i) • In both cases, we arrive at a 
contradiction to our hypothesis. 

Now let us show b). The set S — {gi \ i g {1, . . . , s}, gi ^ (gi, . . . , is a 

system of generators of M , because an element g^ such that gi g (gi, . . . , ^-i) 
is also contained in (gj E S \ 1 < j < i — 1). The fact that this system of 
generators is minimal follows from a). □ 

The following version is an immediate consequence of part a) of the proposi- 
tion. 

Corollary 10 Let N be a graded P -module, let M be a submodule of N, let 
{gi, . . . ,g s } be a minimal homogeneous system of generators of M , and let 
g s +x g iV \M be a homogeneous vector whose degree satisfies the inequality 
deg w (g s+ i) > L ex max{deg H/ (5( i ) | i = 1, . . . , s} . Then {g ly g s+1 } is a min- 
imal system of generators of the module M + (g 8 +i) ■ In particular, we have 
n(M+{g a + x )) = n{M) + \. 



3 Minimal Generators in a Reduced Grobner Basis 

From here on we use the following assumptions. Let if be a field, and let 
P = K\ polynomial ring over K which is positively graded by 

a matrix W g Mat m , n (Z). Then let r > 1, let 5 U . . . , 5 r g Z m , and let M 
be a graded submodule of F — ©[ =1 P(—5i) which is generated by a set of 
non-zero homogeneous vectors {v i, . . . ,v s } . Furthermore, we choose a module 
term ordering a on the monomodule of terms T n (ei, . . . , e r ) in F , and we let 
V = (vi, ...,v s ). 

Our first goal is to describe an algorithm which computes a homogeneous 
cr-Grobner basis of M degree-by- degree and a variant of this algorithm which 
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also yields a minimal system of generators of M contained in V. This part 
is classical and more or less "well-known". Then we make good use of it in 
Theorem 15 for minimalizing reduced Grobner bases. 

To ease the notation, we shall use the following convention: whenever a vec- 
tor gi appears, we write LM CT (<7j) = Cjtje 7 . , where q G K\{0} , where U G T n , 
and where ji G {l,...,r}. For two indices i,j such that ji = 7,-, we let 

_ \cm(U,tj) _ lcm(U,tj) 1 o _ Icmfatj) _ lcmfe.tj) 

** \~>i i/£ 3 3 t 1 3 3 

Theorem 11 (The Homogeneous Buchberger Algorithm) 

In the above situation, consider the following instructions. 

1) Let B = 0, W = V, G = 0, and let s' = 0. 

2) Let d be the smallest degree with respect to Lex of an element of B or 
of W. Form B^ and Wd, and delete their entries from B and W, re- 
spectively. 

3) If B d = 0, continue with step 6). Otherwise, chose a pair G B d and 
remove it from B d . 

4) Compute the S-vector Sy and its normal remainder S'^ = NK^g (Sij) . 
If = 0, continue with step 3). 

5) Increase s' by one, append g s > = 5^- to the tuple Q , and append the set 
{(i, s') I 1 < i < s', 7j = 7 S /} to the set B . Continue with step 3). 

6) If Wd = 0, continue with step 9). Otherwise, choose a vector v G Wd 
and remove it from Wd ■ 

7) Compute v' = NK aj g(v) . If v' = 0, continue with step 6). 

8) Increase s' by one, append g s i = v' to the tuple Q , and append the set 
{(i, s') I 1 < i < s', 7j = 7 S '} to the set B . Continue with step 6). 

9) If B = and W = 0, return the tuple Q and stop. Otherwise, continue 
with step 2). 

This is an algorithm which returns a a -Grobner basis Q of M , where the 
tuple Q consists of homogeneous vectors having non- decreasing multidegrees. 

The proof of this theorem is standard Computer Alg ebra and is for instance 



1 ne proot 01 tin s tneo rem is st andard computer Aig ( 
contained in ( Kreuzer and Robbianol in preparation! ) . 



Remark 12 Let us add some observations about this algorithm. 

a) If we interrupt its execution after some degree d is finished, the tuple Q 
is a <i -truncated Grobner basis of M . Consequently, we can compute 
truncated Grobner bases efficiently. Moreover, in this case it suffices to 
append only the pairs {(i,s') \ 1 < % < s', 74 = j 3 >, deg w (a is >) < Lex ^0} 
to the set B in steps 5) and 8). The reason is that pairs of higher degree 
are never processed anyway, since we stop the computation after finishing 
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degree do. 

b) It is not required that a is a degree compatible module term ordering. 
The reason is that, during the computation of the Grobner basis, only 
comparisons of terms in the support of a homogeneous vector are per- 
formed. Thus these terms have the same degree, and it does not matter 
whether a is degree compatible or not. 

c) The Homogeneous Buchberger Algorithm can also be viewed as a spe- 
cial version of the usual Buchberger Algorithm where we use a suitable 
selection strategy. 

The following variant of the Homogeneous Buchberger Algorithm computes a 
minimal system of generators of M contained in the given set of generators 
while computing a Grobner basis. It provides an efficient method for finding 
minimal systems of generators. 

Corollary 13 (Buchberger Algorithm with Minimalization) 

In the situation of the theorem, consider the following instructions. 

V) Let B = 0, W = V, g = 0, s' = 0, and V min = 0. 

2) Let d be the smallest degree with respect to Lex of an element of B or 
of W. Form B d and Wd, and delete their entries from B and W , re- 
spectively. 

3) If Bd = 0, continue with step 6). Otherwise, chose a pair G Bd and 
remove it from Bd . 

4) Compute the S-vector S^- and its normal remainder S'^ = NR a; g(Sij) . 
If S'^ = 0, continue with step 3). 

5) Increase s' by one, append g s > = S'^ to the tuple Q , and append the set 
{(i, s') | 1 < i < s', 7$ = 7 S /} to the set B . Continue with step 3). 

6) If Wd = 0, continue with step 9). Otherwise, choose a vector v G Wd 
and remove it from Wd ■ 

7) Compute v' = NR ffi g(ti ) . If v' = 0, continue with step 6). 

8') Increase s' by one, append g s > = v' to the tuple Q , append v to the 
tuple Vmm, and append {(i,s') | 1 < i < s', 7« = r y s i} to the set B. 
Continue with step 6). 

9') If B = and W = 0, return the pair (£,V mm ) and stop. Otherwise, 
continue with step 2). 

This is an algorithm which returns a pair (£/, V m i n ) such that Q is a tuple of 
homogeneous vectors which are a o -Grobner basis of M , and V m i n is a sub- 
tuple of V of homogeneous vectors which are a minimal system of generators 
ofM. 

Proof. In view of the theorem, we only have to show that the elements in V m i n 
are a minimal set of generators of M . Since the algorithm is finite, it operates 
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in only finitely many degrees d. Therefore it suffices to prove by induction 
on d that V m i n contains a minimal system of generators of (M<<j) after the 
algorithm has finished working on elements of degree d. 

This is clearly the case at the outset. Suppose it is true for the last degree 
treated before d. Inductively, we can show that the elements of Q continue to 
be contained in the module (M <( j) while we are looping through steps 3), 4), 
and 5) of the algorithm. Namely, every time an element of the form NF CT) g(5'y) 
is added to Q , it is clearly contained in the module generated by the previous 
elements of Q . Furthermore, by part a) of the remark following Theorem 11, 
the elements of the tuple Q form a rf-truncated Grobner basis of (M <£ j) after 
we have finished looping through steps 3), 4), and 5), i.e. when we have treated 
all pairs of degree d. 

Now let Wd = (wi, . . . , we) , and let the numbering of these vectors correspond 
to the order in which they are chosen in step 6). We show that, for each appli- 
cation of steps 6), 7), and 8'), the elements of V m i n continue to be a minimal 
system of generators of the module they generate, and that this module al- 
ways agrees with the one generated by the elements of Q . Furthermore, the 
elements of Q are always a rf-truncated cr-Grobner basis of that module. 

When a new vector v = w% is chosen in step 6), there are two possibilities. 
If v' — in step 7), then v is already contained in the module M' generated 
by the elements of V m j n . Otherwise, the vector v' is not contained in M', since 
the elements of Q are a d-truncat ed cr-Grobner basis and we c an apply the 
Submodule Membership Test (see ( Kreuzer and Robbianol 2000[ ). Proposition 



2.4. 10. a). In that case, the elements of V m i n , together with v, form a minimal 
system of generators of the module M' + (v) = M' + (v') by Corollary 10. 
Moreover, the elements of Q , together with v' , form a <i-truncated cr-Grobner 
basis of M 1 + (v') by Corollary 8. 

Altogether, it follows that, after degree d is finished, the elements of V m i n are 
a minimal system of generators of (M<^) , as we wanted to show. □ 



Remark 14 Let us collect some observations about this algorithm. 

a) If we are only interested in a minimal system of generators of M (and not 
in a Grobner basis), we can stop the algorithm after we have completed 
degree c/ max = max{deg(fj) | 1 < % < s}. In this case it suffices to append 
only the pairs {(i, s') \ 1 < % < s', 7* = 7^, deg w (a is >) < Lex d max } to the 
set B in steps 5) and 8'). 

b) In addition, we could alter step 8') and append the vector v' instead of v 
to the list Vmin- Then V m ; n would still contain a minimal homogeneous 
set of generators of M when the computation ends. These generators 
would not be contained in the initial tuple V anymore, but they would 
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have the additional property that each vector is fully reduced against the 
previous ones. 

The final part of the section is devoted to a result which will be essential for our 
discussion of the minimalization of the critical pairs. Namely, we are going to 
apply the algorithm of Corollary 13 to a reduced Grobner basis and improve it 
significantly in that case. The main differences between both algorithms occur 
in step 7), where it suffices to compare terms instead of computing normal 
remainders, and in step 8), where we append v to both Q and V m i n . 

Theorem 15 (Minimal Generators in a Reduced Grobner Basis) 

In the situation of Theorem 11, let V = (i>i, . . . ,v s ) be the reduced a -Grobner 
basis of M . Consider the following instructions. 

1) Let B = 0, W = V, G = 0, s' = 0, and V min = 0. 

2) Let d be the smallest degree with respect to Lex of an element of B or 
of W . Form Bd and Wd, and delete their entries from B and W , re- 
spectively. 

3) If B d — , continue with step 6). Otherwise, choose a pair G B d 
and remove it from B d . 

4) Compute Slj = NR^g(5j 3 -) . If S'^ = , continue with step 3). 

5) Increase s' by one, append g s i = ■ to the tuple Q , append the following 
set {(i, s') | 1 < % < s', 7, = 7 S /} to B , and continue with step 3). 

6) If Wd = 0, continue with step 9). Otherwise, choose v G Wd and remove 
it from Wd ■ 

7) If LT a (v) = LT a (g) for some g G Q , then replace the element g in Q 
by v. Continue with step 6). 

8) Increase s' by one, append g s > = v to the tuples Q and V m i n; and append 
{(i, s') | 1 < % < s', 7$ = 7 S /} to the set B . Continue with step 6). 

9) If B = and W = 0, return V m i n and stop. Otherwise, continue with 
step 2). 

This is an algorithm which computes a subtuple V m i n of V such that V m i n is 
a minimal system of generators of M . 

Proof. It suffices to show that this procedure has the same effect as running 
the algorithm of Corollary 13 on V. 

First we use induction on d to show that, after we have finished some degree d, 
the tuple Q has the same elements as V< d ■ Every element of Vd is appended 
to Q at some point in step 7) or 8). On the other hand, if an element g s i is put 
into Q in step 5), it has a leading term which is not a multiple of an element 
of V<d. Hence it is swapped out of Q at some point in step 7). 

Next we note that, after we have finished cycling through steps 3), 4), and 5) 
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in degree d, the tuple Q is a d-truncated minimal a-Grobner basis of M<d. 

Now we turn our attention to the loop described in steps 6), 7) and 8). Notice 
that the effect of steps 7) and 8) is independent of the order in which we 
choose the elements v G Wd in step 6). Hence we can assume for the purposes 
of this proof that we always choose the vector v in Wd which has the minimal 
leading term with respect to a. With this assumption, we show inductively 
that when we run steps 7) and 8) for some element v G Wd , at each point the 
elements in Q are a minimal cr-Grobner basis of the module they generate, 
and the elements of V m i n are a minimal system of generators of that module. 

For the induction step, we have to consider two cases: either v is swapped 
into Q in step 7) or appended to both Q and V m i n in step 8). In the first case, 
it suffices to show that the module generated by the elements of Q does not 
change when we perform the swap, i.e. that the difference v — g is contained 
in this module. This follows from the observations that LT a (v — g) < a LT a (v) 
and all elements v in V such that LT a (v) < a LT a (v) are already in Q . Since 
v — g — > , we have v — g — > . In the second case, it is clear that Q continues 
to be a minimal Grobner basis of the module it generates by Corollary 8, 
and V m in continues to be a minimal system of generators of that module by 
Corollary 10. 

Finally, we note that in step 8) we can append v to Q without passing to the 
normal remainder, since v is an element of a reduced Grobner basis and thus 
irreducible. □ 



Remark 16 Let us make some observations about the preceding algorithm. 

a) The proof of the proposition shows that the algorithm reconstructs the 
given reduced Grobner basis inside Q, and that Q<d has the same ele- 
ments as V<d after some degree d is finished. 

b) Moreover, we note that in step 4) it is not necessary to compute the 
normal remainder NR ff) g(5jj) . Rather, it suffices to perform a full leading 
term reduction. 

c) The different elements NR^g(5i 3 -) computed in step 4) and the elements 
v G Vd which are swapped into Q by step 7) are in 1 — 1 correspondence, 
since every new element computed in step 4) must have a new leading 
term in the leading term module of M . This new leading term must be 
the leading term of an element in the reduced Grobner basis, hence it is 
swapped. 



12 



4 Minimalizing the Critical Syzygies 



In this section we continue to use the assumptions and notation of the previous 
section. If we look at Theorem 11 and its proof, we can see that instead 
of treating all pairs (i,j) such that is contained in the set of critical 
syzygies £ , it would be enough to treat those pairs corresponding to a subset 
6 C £ which is a minimal system of generators of Syz P (citie 7l , . . . , c s t s e 7 J . 

In order to find G, we obse rve that the applicat i on of two of the rules for 
killing critical pairs given in (|Oebauer and Molleil . fl987h produces a minimal 



Grobner basis of the module Syz P (citie 71 , . . . , c s t s e ls ) contained in the set S. 
From this we derive the idea to find by applying Theorem 15. We need the 
following definition. 



Definition 17 On the set of terms T n {e\, . . . ,e s ) in ©| =1 P(—di) we define 
a relation r by letting 



t Ei > T t' Sj ^ | 



t ti 6.-^ ^> & t tj , Or 

t U e 7i = t' tj e 7 . and i > j 



for t,t' G T n and i,j € {!,...,«} . As in ( Kreuzer and Robbianol 2000) 



Lemma 3.1.2, it follows that r is a module term ordering. It is called the term 
ordering induced by the tuple (tie 71 , . . . , t s e ls ) and by a. 



By ( Kreuzer and Robbianol EoOOj ) . Proposition 3.1.3, the set S is a r-Grobner 



basis of the module Syz P (citie 7l , . . . , c s t s e la ) . Moreover, cr^ is a homogeneous 
element of ffi? =1 P(— dj) of degree deg l y((7y) = deg(lcm(tj, tj)) + S yi . For all 
i, ?' 6 {1, . . . , s} , we let Uj = lcm ( f " f j) . n ow the main result of lGebauer and Mollei 
reads as follows. 



Proposition 18 Consider the following instructions. 

Rule 1 . Delete in £ all elements Oj^ such that there exists an index i in the 
set {1, . . . , j — 1} such that t^ divides t^j . Call the resulting set £'. 

Rule 2. Delete in £' all elements such that there exists an index j in the 
set + . . . ,k — 1} such that t^j properly divides . Call the resulting 
set £" . 

Rule 3. Delete in S" all elements o~ij such that there exists an index k in 
the set G {j + 1, . . . , s} such that t ik properly divides tij and tj k properly 
divides tji . Call the resulting set . 

Then the set still generates Syz P (citie 7l , . . . , c s t s e ls ) . 
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Remark 19 Let us interpret the previous proposition in another way. For 
1 < i < j < s such that 7, = jj, we have LT r ((Ty) = tjj£j. Hence Rules 1 
and 2 can be restated as follows. 



Rule 1'. Delete in E all elements such that there exists an element a^j 

such that LT T (cry) is a proper multiple of LT T (oVj). 
Rule 2'. If, among the remaining elements, there are elements a^, a^j such 

that UT T ((Jij) = LT r (crj'j), then delete the one having the larger index 

max{i, i'} . 

From Rules 1' and 2' it follows that the set E" is a minimal r-Grobner basis 
of the module Syz P (citie 7l , . . . , c s t s e 7s ), i.e. the leading terms of the elements 
of E" minimally generate the leading term module. 



In general, it is not true that E" is a minimal system of generators of the 

module Syzp fcjtie^ c.,t.,e, v „'), as ou r next example shows. (For another 

example, see ([Gebauer and Moiled . 119871 ). 3.6.) 



Example 20 Let P = Q[x,y,z] be standard graded, let r = 1, s = 4 and 
t\ = x 3 z 2 , t 2 = x 3 y A , £3 = y 5 z 2 , £4 = x 2 y b z. Then we get a 12 = y A e\ — z 2 E2, 
&13 = y 5 e 1 -x 3 e 3 , a u = y 5 e 1 -xze 4 , 023 = yz 2 e 2 -x 3 e 3 , o 2 a = yze 2 -xe 4 , and 
C34 = x 2 Es — ze<±. By applying Rules 1 and 2, we get the minimal r-Grobner 
basis E" = {0-12,^4,^34,^13} of Syz P (ti, t 2 , t 3 , t 4 ) , since LT T (cr 23 ) = LT r (cri 3 ) 
and LT T (o"i 4 ) = z ■ LT r (cx 2 4). Now we use Rule 3 and find E'" = E", but E w 
is not a minimal system of generators of Syz P (£i, t 2 , t 3 , t^j , since we have 



Before continuing, let us introduce a new notion. If we have an element err- 
and perform a reduction step — *-l dt'si + c"t"e^ , where c, c', c" G K and 
t,t',t" G T n , we call this a head reduction step. (Notice that the j'-indices 
have to match!) Similarly, we can define a tail reduction step as follows: 

<Tij c -^-i c't'si' + c"t"ej . It is clear that a tail reduction step does not change the 
leading term of the element. 

Proposition 21 The set E = {— Cj ■ Oij \ G E"} is the reduced r-Grobner 
basis of the module Syz P (citie 7l , . . . , c s t s e ls ) . 

Proof. Since passing from E" to E is equivalent to normalizing the leading 
coefficients, and since E" is a minimal r-Grobner basis, it remains to show that 
no tail reductions are possible among the elements of E . But if we perform a 
tail reduction on some element of E , we get an element of the form cto~i'j such 
that i' <%. Here we have to have t = 1, since (Tjj is part of a minimal Grobner 
basis. Now we obtain a contradiction to the minimality of i in Rule 2'. □ 
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Remark 22 Let us apply the algorithm of Theorem 15 to the set E. We make 
the following observations. 

a) A pair of pairs, i.e. a critical pair between two elements a^, <7,y yields an 
S- vector ^((ij),^',^)) = c ^ a ij ~ dt' ovy such that c,c' & K and t,t' G T n 
and j = f , since the two leading terms have to cancel. Without loss 
of generality, let i < i' . Then the result is ctau' for some c E K and 
t G T n . The degree of such a pair of pairs is 

deg w (S {{itjWd)) ) = deg w (t) + deg PF (a ii ,) 

= deg w (^^)+deg w ( £j ) 
= deg w (lcm(t i ,t i /,t j )) + 5 7j 

b) During the course of the algorithm, a new Grobner basis element can 
only be obtained from a pair of pairs if t — 1 . This is equivalent to 
gcd^ti/j) = 1. 

Now we are ready to optimize the minimalization of the critical syzygies. To 
ease the notation, we shall minimalize the set E" instead of E. The lack of 
the normalization of the leading coefficients is clearly of no consequence. We 
need the following lemma. 

Lemma 23 Let \ <i<j<m<s and i' G {1, . . . ,j — 1} \ {i} . Suppose 
there are terms t,t',t" G T n \ {1} such that = + t a^j = t' a im — t" o~j m 
and o~i' m = tai'j +t" Uj m . Then t, t! , and t" are pairwise coprime. 

More precisely, given k G {1, . . . , n} , we define a = deg 2 , K (t i ) ; a' = deg^^), 
(3 = deg XK (tj) , and 7 = deg XK (t m ) . Then one of the following four cases 
occurs. 

1 ) We have a = 7 > (5 and a > a' . 

2) We have a' = (5 > 7 and a' > a. 

3) We have a = a' > (3 and a > 7 . 

4 ) We have a = a' = f3 > 7 or a = /? = 7 > a' or a' = = 7 > a. 

Proof. Comparing coefficients in the given equations yields the following equal- 
ities \cm(ti,tj) = lcm(tj,tj/) = lcm(tj/,t m ) = t lcm(tj/,t,-) = t' lcm(tj,t m ) = 
t" \cm(tj,t m ) . Thus the exponent of x K in these terms satisfies max{o;,/3} = 
max{«, a'} = max{a', 7} = deg XK (t) + maxja', (3} = deg^t') + maxja, 7} = 
deg XK (t") + max{/9,7}. We distinguish the following four cases. 

Case 1: Suppose that x K divides t. In this case, max{a,a'} > m&x{a',(3} 
yields a > a' and a > (3. Then a = max{a, a'} = max{«', 7} shows a = 7, 
i.e. we have the inequalities stated in case 1) of the claim. Furthermore, it 
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follows that 7 = max{a,7} = max{/9,7}, i.e. that x K divides neither t' 
nor t" . 

Case 2: Suppose that x K divides t! . In this case, max{a, a'} > max{a:,7} 
yields a' > a and a' > 7. Then maxjo;,/?} = max{«, «'} shows a' — (3, i.e. 
we have the inequalities stated in case 2) of the claim. Furthermore, it follows 
that (3 = max{a',P} = max{/3, 7}, i.e. that x K divides neither t nor t" . 

Case 3: If x K divides t" , we argue analogously and obtain the inequalities 
stated in 3) as well as the fact that x K divides neither t nor t' . 

Case 4: If x K divides neither t nor t! nor t" , an easy case-by-case argument 
yields the possibilities listed in 4). □ 

Proposition 24 (Minimalization of the Critical Syzygies) 

Let E" be the r-Grobner basis of Syz P (ciiie 7l , . . . , c s t s e ls ) defined in Propo- 
sition 18. Consider the following instructions. 

1) Let B* = 0, W = E" ; A = ®, and = 0. 

2) For all 01^0^ j G S" such that l<i<i'<j<s, form the S-vector 

= t On 1 , where t G T n . If t — 1 , append aa> to B* . 

3) Let d be the smallest degree with respect to Lex of an element of B* or W. 
Form B* d and Wd, and delete their entries from B* and W , respectively. 

4) If B* d = , continue with step 11). Otherwise, choose an element G B* d 
and remove it from B* d . 

5) If LT T ((jjj) G LT T (Ad) , then continue with step 4)- 

6) If LT T (o"ij) = LT r (o"j'j) for some element a^j G Wd, then remove a^j 
from Wd, append it to A, and continue with step 4)- 

7 ) Find a^j G A<d such that t^ is a multiple of tjy . Then perform the 

head reduction step Oij —^iou, where i 6 T", where k = min{i,i'}, 
and where £ = max{i, i'} . If t 7^ 1 , continue with step 4 )■ 

8) If LT T (ake) £ LT T (^4^) ; then continue with step 4)- 

9) If LT T (<7(^) = LT T (c>V£) for some element oyt G Wd, then remove the 
element a^e from Wd, append it to A, and continue with step 4)- 

10) If uu G B* d , then delete au in B d and continue with step 7), applied to 
this element. Otherwise continue with step 4)- 

11) Append Wd to A and to ©. 

12) If B* = and W = 0, return 6 and stop. Otherwise, continue with 
step 3). 

This is an algorithm which computes a subset 6 C E" such that is a 
minimal system of generators of Syz P (citie 71 , . . . , c s t s e ls ) . 

Proof. It suffices to show that the given instructions define an optimization of 
the application of Theorem 15 to the set E". The tuple A corresponds to Q 
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there, 6 corresponds to V m in, and B* corresponds to B . 

The first significant difference occurs in step 2). Instead of producing the pairs 
of pairs inductively each time we find a new Grobner basis element, we pre- 
compute them all at once. This is possible, since we know from Theorem 15 
that we are merely recomputing the Grobner basis £". Moreover, we do not 
store the pairs of pairs, but the ^-vectors they generate, and we do not store 
^-vectors which are clearly useless by part b) of the remark following Propo- 
sition 21. 

The main difference occurs in steps 5) through 10). Instead of computing the 
normal remainder of the S-vector, we perform leading term reductions only 
and check the result after each reduction step. When we choose an element 
in step 4), it is not contained in Ad, since if an element au is appended to A 
in step 11) and cannot be contained in E* d by step 10). But the element 
could have a leading term in UI a (Ad) without being contained in Ad- We 
claim that, in this case, we know <7jj , i.e. that produces no new 
Grobner basis element. 

To prove this claim, we first note that clearly A is a subtuple of W at all 
times. Since the elements of W are fully interreduced, the tail of aij cannot be 
a leading term of an element of Ad - On the other hand, if LT T (o"jj) = LT T (oVj) 
for crj'j G Ad , then the leading term of the result of the reduction of by 0Vj 
is the tail of . Hence can be tail reduced using A<d ■ By applying the 
same argument to the result of this tail reduction step, we conclude that after 
several tail reductions using A<d, we reach an element of Ad, and the claim 
follows. 

The next possibility for is that it is head irreducible with respect to A. 
In this case its leading term is equal to LT T (<7j/j) for some avj G Wd- Now 
Theorem 15 says that we should put NR,-^^) into A and later swap it 
for <7j/j. But, as we just saw, we can tail reduce using A<d until we 
reach a^j . Thus the normal remainder is a^j and is put into A immediately, 
i.e. without actually performing the tail reductions and without a later swap. 

The last possibility for LT T (<7jj) is that it can be reduced using A<d- This 
reduction step is performed in step 7). Let us discuss the possible outcomes. 

If the result is of the form ta k £ with t G T™ \ {1}, then au has a lower 
degree and satisfies au — ► , because A contains a truncated Grobner basis. 
Consequently, we have and step 4) of 15 tells us to try the next 

S-vector. 

If the result of the head reduction step has one of the new leading terms 
provided by the elements of Wd, we notice this in step 8) or 9). In the first 
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case, the element of Vd has already been swapped into A and nothing needs 
to be done. In the second case, we perform the swap in step 9). 

If the result is an element au of degree d which can be further head reduced, 
we check in step 10) whether au G B* d . In that case and au have the same 
reductions and it suffices to treat <7ke in step 7). Otherwise, we claim that (?ki 
is one of the elements of B d which has been dealt with already, i.e. that we 
can go back to step 4) and treat the next element of B* d . 

To prove this claim, we first use aij G B* d in order to write = t' a im + 1" <Jj m 
with t',t" G T n \ {1} and j < m < s . Secondly, by step 7), we have the 
equality = ta^j ± a ki) where Gu = ±oyj and t G T n \ {1}. By looking at 
the coefficient of ej in the equation = t' a irn — ta^j — t" <Jj m , we see that 
t \cm(ti',tj) = t" lcm(tj, t m ) . This term is a multiple of and of t m . Hence it 
is of the form t lcm(tj/, t m ) for some t G T n , and we have = t' o im — t o"j/ m . 
If t 7^ 1, then ou is a pair of pairs, i.e. it is either in B d or it is one of the 
elements of B d treated before. Hence the claim follows if we can show that 
t — 1 does not happen. 

Suppose that t — 1 . Then we are in the situation of the lemma. Since the 
conditions of steps 8) and 9) did not apply, it follows that <Jkt can be further 
head reduced using A <( i- Hence there exist u,u' G T n and j' < m&x{i,i'} 
such that cTj/j = u a^j* + u' oyi and u ^ 1 or u' ^ 1 , depending on whether 
i > i' or i < i' . 

Now we show that u' ^ 1 is impossible. We use the notation of the lemma 
and let 5 = deg XK (tj>), where x K is one of the indeterminates occurring 
in t, i.e. where case 1) of the lemma holds. Then the equation lcm(tj/,tj) = 
u lcm(tj/, tj>) = u' \cm(ti,tji) shows m&x{a,a'} > m&x{a,5}. This implies 
a' > a and a' > 8, in contradiction to case 1) of the lemma. Similarly, we 
can show that u ^ 1 is impossible. This concludes the proof of the claim. 

Altogether, it follows that steps 5) - 10) implement the full reduction of al- 
together with the swapping procedure of step 7) of 15. Hence the remaining 
elements of Wd are precisely the minimal generators of degree d we are looking 
for, and they have to be appended to 6 in step 11). □ 

Let us apply this algorithm in the situation of Example 20. 

Example 25 Our task is to minimalize W = S" = {a± 2 , o"i3, cr 2 4, o" 34 }, where 
we have deg w (a 12 ) = 9, deg w (cri 3 ) = 10, and deg w ((T24) = deg w (a 34 ) = 9. 

In step 2), the algorithm constructs the set B* . The pair of pairs ((2, 4), (3, 4)) 
yields 5 , (( 2 ,4) 1 (3 1 4)) = za 2 i — xcr 34 = —yz 1 e 2 + x 3 £3 = 023, and this is the only 
element of B* . Notice that it has degree 10. 
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In step 3), the algorithm starts to operate in degree d — 9. Since Bg — 0, it 
appends a 12 , a 24 , and 034 to „4 and 6 in step 11). 

Next we process degree 10. In step 4), we choose <7 2 3 G B\ and set B\ = 0. 
Then, in step 6), we find LT r (<j 2 3) = x 3 e 3 = LT T (o" 13 ), where a± 3 G Wio- 
Hence a± 3 is removed from Wio and appended to A in step 6). 

Thus we have B* = and W = at this point, and step 12) returns the set 
© = {C12, o"24, 034} • We note that this is the correct answer, and there is an 
improvement over the application of Proposition 18 coming from the fact that 
in step 6) we merely check LT T ((7jj) G LT r (Wd) rather than Uij G Wd- 

The following example provides a case where it is actually necessary to do 
one head reduction step in 7) in order to find a previously undiscovered non- 
minimal critical syzygy. 

Example 26 Let P = Q[xi, . . . , x 5 ] be standard graded, let r = 1 and s = 4. 
The terms t\ = x\x\x 4 x\ , t 2 = xfx 2 x 4 xf , t 3 = x^x^Xg , and t 4 = 2^X3X5 yield 
the critical syzygies o\ 2 = x\x1e\— x 2 x 3 e 2 , (T13 = xfei— X4Xg£3, cr u = xfx|£i — 

X2X4&4 , (T23 = X2X3^2 _ ^4^5^3 , 0"24 = ^3^2 ~ X2X4^4 , and (T34 = xf^ — X2^4 . 

Here steps 1) and 2) of Proposition 18 discard a 23 and <t 14 , because we have 
LT T (<j 2 3) = x 4 X5£ 3 = x\ LT T (a- 13 ) and LT r (<r 14 ) = X2X 4 £4 = x 2 LT T (o- 24 ) . Thus 
we have S" = {a i2 , a 13 , <t 24 , a 34 } . We note that we have deg w (a 12 ) = 21, 
deg H /(o"i 3 ) = 19, and deg w (<j 24 ) = deg w (au) = 20. But E" is not minimal, 
since we have o\ 2 = x\a\ 3 — x 2 a 24 + X4O34. 

Now we apply our algorithm. In step 2), we have to compute 5 , (( 2 ,4) 1 (3 1 4)) = 
x 2 o"24 — = ^22^3^2 — ^4^5^3 = °"23 • Thus <t 2 3 is appended to B* . It has 

degree deg w (cr 23 ) =21. No further pairs of pairs are found. 

In step 3), the algorithm starts to operate in degree d — 19. We have B\ 9 = 
and W19 = (cr 13 ). Thus we append a 13 to A and 6 in step 11). Next we pass 
to degree d — 20. We still have B 20 = 0, but now we get W20 = (024, 034) • In 
step 11), 024 and 034 are put into A and 6. 

When we start processing degree d = 21 , we have to choose a 2 3 G B 21 and set 
B 21 = in step 4). The leading term LT r (a 2 3) = x 4 X5£ 3 is not equal to one of 
the leading terms of the elements of A 2 i or W21 • But we can perform a head 
reduction step in 7), namely C23— o"i2- Here step 8) does not apply, but 
in step 9) we have LT r (o"i 2 ) G LT r (W2i). Thus we continue by removing o\ 2 
from W21 and appending it to A. 

Finally, we get B* = and W = 0. The algorithm returns = {013, <t 2 4, 034} . 
As mentioned above, the non-minimal critical syzygy o\ 2 was discovered after 
one head reduction step in 7). 
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5 An Optimized Buchberger Algorithm 

In this section we combine the results obtained so far. We continue to use 
the notation and conventions of the previous sections. In particular, we let 
P = K[xi, . . . ,x n ] be a polynomial ring over a field K which is positively 
graded by a matrix W G Mat m>n (Z), and we let M be a graded sub module 
of a graded free P-module F = (& r i=1 P{— 8A which is generated by a tuple 
V — (vi, . . . ,v s ) of homogeneous vectors. Furthermore, we let a be a module 
term ordering on T n (ei, . . . , e r ) . 

In the following theorem the sets of critical pairs corresponding to the sets 
of critical syzygies considered earlier are denoted by the normal letters corre- 
sponding to their calligraphic versions. 

Theorem 27 (Optimized Buchberger Algorithm) 

In the above situation, consider the following sequence of instructions. 

1) Let W = V, A = ; B = 0, B* = 0, Q = 0, and let s' = 0. 

2) Let d be the smallest degree w.r.t. Lex of an element of B or W . Form 
B d , B* d , W d , and delete their entries from B , B* , and W, respectively. 

3) Apply MinPairs(A,S d ,S^). 

4) If Bd = 0, then continue with step 7). Otherwise, choose a pair 
in Bd , delete it from Bd , and append it to A. 

5) Compute Sij and S[j = NR at g(Sij) . If S'^ = 0, then continue with 4)- 

6) Increase s' by one, append g s > = S'^ to Q, perform Update(£>, B*, g s >) , 
and continue with step 4)- 

7) If Wd = then continue with 10). Otherwise, choose v G Wd and delete 
it in W d ■ 

8) Compute v' = NR CT) g(t>) . If v' = 0, continue with step 7). 

9) Increase s' by one, append g s > — v' to Q and perform Update(5, B* , g s i) . 
Then continue with step 7). 

10) If B = and W = 0, then return Q and stop. Otherwise, continue with 
step 2). 

Here the procedure Update(i?, B*,g s /) is defined as follows. 
Ul) Form the set C = {(i, s') | 1 < % < s', 7, = 7s'}- 

U2) Delete from C all pairs (j, s') for which there exists an index i in the set 

{1, ... ,j — 1} such that t s 'i divides t s 'j . 
US) Delete from C all pairs (i, s') for which there exists an index j in the 

set {i + 1, . . . , s' — 1} such that t s 'j properly divides t s n . 
V 4) Find in C all pairs (i,s') and (j,s') such that 1 < % < j < s' and 

such that gcd(ti S i,tj S >) = 1. For each of these, check if is already 

contained in B* and append it if necessary. 
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U5) Append the elements of C to B and stop. 

Furthermore, the procedure MinPairs(A, B d , B d ) is defined as follows. 

Ml) If B^ — ; then stop. Otherwise, choose a pair in B^ and remove 

it from B* d . 

M2) If tji = tji' for some pair G A, then continue with step Ml). 

MS) If t^ = tj^ for some pair G B d , then remove this pair from B d and 

append it to A. Continue with step Ml). 
M4) Find G A such that tji> divides tji. Let k = min{i,i'}, and let 

£ = max{i,i'}. If gcd(tij,ti'j) ^ 1, then continue with Ml). 
M5) If tgk = tik' for some pair (k',£) G A, then continue with Ml). 
M6) If = tik> for some pair (k', £) G B d , then delete this pair in B d , append 

it to A, and continue with Ml). 
Ml) If (k,£) G B d , then delete (k,£) in B d and continue with M4), applied 

to this pair. 
M8) Continue with step Ml). 

Altogether, we obtain an algorithm which computes a tuple Q whose elements 
form a homogeneous a -Grobner basis of M . Moreover, the set of pairs which 
are treated at some time in steps 4) - 6) of the algorithm corresponds to a 
minimal system of generators of the module Syz P (c 1 t 1 e 7l , . . . , c s it s ie ls ,) . 

Proof. The main algorithm of this theorem agrees with the Homogeneous 
Buchberger Algorithm (see Theorem 11), except for the introduction of the 
procedure MinPairs(A, B d , B d ) in step 3) and the alteration of the enlarge- 
ment of B in steps 5) and 8) of Theorem 11 which is now performed by the 
procedure Update(i?, B*, g s >) . 

The foundation for these changes is the material presented above, especially 
Proposition 24. In steps 4) - 6) we want to treat only those pairs (i,j) for 
which the corresponding elements are contained in the minimal system of 
generators 9 of the graded P-module Syz P (citie 7l , . . . , c s it s ie ls ,) . 

Procedure Update(i?, B*, g s >) applies Rules 1) and 2) of Gebauer-Moller in 
steps U2) and U3), respectively. Moreover, notice that step U4) computes all 
pairs of pairs which satisfy the condition of part b) of Remark 22, and stores 
the pairs corresponding to the resulting S- vectors in B* . 

Thus, in order to minimalize the critical pairs we process, we need to apply 
Proposition 24 to the set of critical syzygies corresponding to the set of critical 
pairs B , where we can refrain from computing the pairs of pairs, because they 
have already been generated and stored in B* . This task is performed by the 
procedure MinPairs(A, B d , B d ) . Its steps Ml) - M8) are easy translations of 
steps 4) - 10) of Proposition 24 into the language of pairs. Notice that we 
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have LT T (a"ij) = LT T (o"^) if and only if j — £ and tji = t^- Altogether, 
Update(S, B*, g s i) and MinPairs(A, Bd, B*j) make sure that only the pairs 
corresponding to are treated at some point in steps 3) - 6). 

Finally, we remark that A is used to keep track of the pairs for which <7jj 
is in that part of the minimal r-Grobner basis E" of Syz P (citie 7l , . . . , c s it s ie ls ,) 
which has been computed so far. Thus it is updated when a non-minimal ele- 
ment of E" is found in step M3) or step M6), and when a pair corresponding 
to an element of is chosen for treatment in step 3). □ 

Let us illustrate the performance of this algorithm by a simple example. It 
shows that cases like Example 20 occur naturally during actual Grobner basis 
computations. 

Example 28 Let P = Q[x,y,z] be standard graded, let a = DegLex, let 
r = 1 , and let M C P be the homogeneous ideal generated by the polynomials 
v 1 = x 3 z 2 + x 2 y 2 z, v 2 = x 3 y s , and v 3 = y 10 z 2 . Then the leading terms are 
ti = x 3 z 2 , t 2 = x 3 y 8 , and t 3 = y 10 z 2 . Let us follow the steps of the Optimized 
Buchberger Algorithm. 

The first degree is d — 5 . Since B 5 = , the first actions are to choose v\ G W5 
in step 7) and append g 1 — V\ to Q in step 9). Then we continue with d — 11 
and choose v 2 G Wn in step 7). Since v' = NR . i g(t> 2 ) = v 2 , we append^ = v 2 
to Q in step 9) and update the set of pairs. The result is B = {(1,2)} and 
B* — 0. Now we have to treat the degree d = 12. Notice that the degree 
of the pair (1,2) is 13. Hence B\ 2 = and we have to choose t> 3 G W± 2 in 
step 7). Since v' = NR ff] g(t) 3 ) = t> 3 , we append g 3 = t> 3 to Q in step 9) and 
update the set of pairs. In step Ul), we form C = {(1, 3), (2, 3)}. In step U2), 
we obtain t 31 = x 3 = t 32 , and therefore (2,3) is deleted in C. The result is 
B = {(1,2), (1,3)} and B* = 0. This completes degree 12, and we continue 
with degree 13. 

We choose the pair (1,2) in step 4) and append it to A. Then we compute 
S 12 = y 8 gi — z 2 g 2 = x 2 y w z and S[ 2 = NR - i g(Si 2 ) = x 2 y w z. Thus we have a 
new Grobner basis element g 4 = x 2 y 10 z and need to update the pairs again. 
In step Ul), we form C = {(1,4), (2,4), (3,4)}. Step U2) does not apply, but 
in step U3) we remove the pair (1,4) from C, since t^ 2 = x properly divides 
£41 = xz. Now we check that t 2 ± = y 2 z properly divides t 23 = y 2 z 2 and 
t 3 4 = x 2 properly divides t 32 = x 3 . Hence the pair (2, 3) is appended to B* . 

At this point we have finished degree 13, and we have the following situation: 
A = {(1,2)}, B = {(1,3), (2,4), (3,4)}, B* = {(2, 3)}, Q = { 9l , . . . , g 4 }, and 
s' — 4. The next degree is d — 14, where we have to deal with the pairs in 
514 = {(2,4), (3,4)}. Since B* 4 = 0, we choose (2,4) in step 4) and append 
it to A. Then we compute S 24 = and continue by choosing (3,4) in B 14 
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and adding it to A. Again S34 = , and degree 14 is finished. 

Now we start degree 15 by performing MinPairs(A, Bi$, B^) , where we have 
A = {(1,2), (2,4), (3,4)}, B 15 = {(1,3)} and B* 5 = {(2,3)}. In step Ml), we 
choose (2,3). In step M3), we discover t 32 = x 3 = i 3i , where (3,1) G B 15 . 
Hence (1,3) is deleted in B\<~, and appended to A. Then the procedure is 
finished, and the facts that B^ = as well as W15 = allow us to return Q 
and stop. 

As in Example 20, we have found one useless pair, namely the pair (1, 3) 
in degree 15, which would not have been discovered by the Gebauer-Moller 
Installation, and which we were able to discard by a simple combinatorial 
check. 



Remark 29 Let us discuss the efficiency of the algorithm of Theorem 27. 

a) Steps U2) and U3) of this algorithm correspond to Rules 1) and 2) of 
the Gebauer-Moller installation. However, Rule 3 is not performed by the 
procedure Update(. . .) , but by step M2) of the procedure MinPairs(. . .) . 
In fact, step M2) gets rid of more pairs than Rule 3, because Rule 3 
requires (i, j) G B d D B d , whereas we only need a pair G B d such 
that LT T (<7jj) = LT T ((7j/j) for some G B d . 

b) A potential drawback of our approach is that the number of pairs of 
pairs considered in step U4) is quadratic in the number of elements of C 
surviving steps U2) and U3). But that number is usually fairly small. 
Hence the cost of U4) and the cardinality of B* tend to be rather small. 
On the other hand, we do not need to check Rule 3 for all elements of 
the list B which is usually rather long. Our experiments suggest that, on 
average, the overhead of the two approaches is comparable. 

c) Our procedure MinPairs(. . .) is very efficient in treating the elements 
of B* d . Each time we loop through steps M2) - M8), we delete one pair 
in B d , and B d is never enlarged. In practice, we find that the lists B d 
are generally small. Hence our algorithm harnesses the full power and 
efficiency of the Gebauer-Moller installation, while it simultaneously kills 
all unnecessary pairs at a comparatively small cost. 



6 Experimental Data and Conclusions 

In this section we want to provide the reader with some experimental nu- 
merical data which illustrate the performance of the Optimized Buchberger 
Algorithm 27 as well as technical observations coming from an implementation 
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in an experimental version of the "CbCbA 5" library in C++. 

In the following table, we compare the application of Rules 1) - 3) of Proposi- 
tion 18 to our procedures Update(. . .) and Minpairs(. . .) in Theorem 27, i.e. 
to the algorithm of Proposition 24. Let us point out that our procedure al- 
ways minimalizes the critical pairs, independent of the order of the underlying 
terms. (Non- minimal critical pairs are recognized at different steps, though.) 
For the Gebauer-Moller installation, however, the number of undiscovered non- 
minimal critical pairs depends strongly on this order. 

To aid the reader in understanding this table, let us explain the meaning of 
the symbols. 

• is the cardinality of the reduced Grobner basis of the corresponding 
ideal. 

• is the total number of pairs, i.e. = (^f^) ■ 

• is the number of pairs surviving Rules 1) and 2), i.e. the cardinality 
of the reduced Grobner basis of pairs. 

• B is the number of pairs killed by Rule 3), the Gebauer-Moller "Back- 
wards" criterion. 

• M23 is the number of pairs killed by steps M2) and M3) in Theorem 27. 

• M48 is the number of pairs killed by steps M4) - M8) in Theorem 27. 

• Gain = M23+M48— B , i.e. the number of newly discovered non-minimal 
critical pairs. 

• #(©) is the cardinality of a minimal system of generators of the syzygies 
of the leading terms. Hence we have #(6) = #(£") - M23 - M48 





#(G) 


#(£) 


#(£") 


B 


M23 


M48 


Gain 


#(6) 


T~51 


83 


3,403 


250 


7 


7 








243 


Twomat3 


109 


5,886 


741 


15 


26 


1 


12 


714 


Alex3 


211 


22,155 


684 


54 


56 


1 


3 


627 


Gaukwa4 


267 


35,511 


1,772 


101 


113 


3 


15 


1,656 


Kinl 


306 


46,665 


3,411 


70 


172 





102 


3,239 


Wang (Lex) 


317 


50,086 


1,457 


60 


61 


7 


8 


1,389 


Cyclic 7 


443 


97,903 


2,651 


17 


17 








681 


Hairer-2 


506 


127,765 


5,305 


150 


152 


4 


6 


5,149 


Hom-Gonnet 


854 


364,231 


11,763 


587 


648 


27 


88 


11,088 


Mora-9 


4,131 


8,530,515 


46,395 


1,930 


1,914 


23 


7 


44,458 



The rows of this table correspond to standard examples of Grobner basis 
computations. A file containing a description of every example can be down- 
loaded at 
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ftp: //co coa . dima . unige . it/papers/ CaboaraKreuzerRobbiano03 . cocoa 



Moreover, a file containing the list of leading terms of the reduced Grobner 
basis for each example can be downloaded at 

ftp : //cocoa . dima . unige . it/papers/CaboaraKreuzerRobbiano03_2 . cocoa 

Technical note: In the well-known example "Cyclic 7" we have homogenized 
using a new smallest indeterminate (see the file mentioned above). 

For the reader who would like to run his own tests, we note that , , 

and #(©) are invariants of the reduced Grobner basis. But the effect of both 
the Gebauer-Moller installation and our Optimized Buchberger Algorithm de- 
pend strongly on the order in which the elements of £ are produced during 
a Grobner basis computation. For instance, this means that it depends on 
the chosen selection strategy. In our implementation pairs are kept ordered in 
increasing DegLex ordering, reductors are kept in the same order they are pro- 
duced, reductors of the same degree are kept interreduced, and the reduction 
strategy is full reduction. 

The following table shows some timings. It compares Singular 2.0.0 with the 
current experimental version of CbCbA 5 using the GM and CKR pair han- 
dling algorithms. Timings are in seconds for Linux running on an Athlon 
2000+ CPU with 1.5GB RAM. All computations are over the rationals where 
the timings of the base field operations in Singular and CoCoA seem to be 
comparable. 

Technical note: The reason why we include a comparison with Singular is an 
explicit request made by a referee, who suggested comparing our timings with 
"another efficient implementation" . The table below indicates that both Sin- 
gular and CoCoA 5 have efficient implementations of the Buchberger algorithm, 
and that our new algorithm has at least the same efficiency. 





Singular 2.0.0 


CbCbA5 GM 


GoCbA5 CKR 


T"51 (Lex) 


149.32 


7.28 


7.14 


Twomat3 


1.21 


8.66 


8.50 


Alex3 


«1 


0.54 


0.56 


Gaukwa4 


80.30 


99.31 


98.57 


Kinl 


407.09 


89.25 


87.41 


Wang (Lex) 


>1200 


382.86 


379.31 


Cyclic 7 


>1200 


76.61 


76.65 


Hairer-2 


79.36 


141.83 


139.76 


Hom-Gonnet 


3.97 


4.55 


4.95 


Mora-9 


30.53 


86.17 


89.75 
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Conclusions 



First of all, let us collect some technical observations based on our implemen- 
tation of the Optimized Buchberger Algorithm. 

a) When we apply Rules 1 and 2 of Proposition 18, the remaining set of 
pairs E" is usually almost a minimal system of generators of the mod- 
ule Syz P (citie 7l , . . . , c s it s ie ls ,) . Thus both Rule 3 and our algorithm kill 
comparatively few pairs. Nonetheless, over the rationals (or other costly 
fields), the saving is worthwhile because the treatment of each single pair 
can take a long time. 

b) Steps M5) - M7) in the Optimized Buchberger Algorithm are indepen- 
dent. Hence it is possible to order them in such a way that the com- 
putational cost is minimized. This may be important if there is a large 
number of elements in B^ to be processed, since the operations may have 
substantially different computational costs. 

c) All operations in our procedures Update ( . . . ) and MinPairs( . . . ) have 
been greatly eased by memorizing the terms tij, tji and lcm(tj, tj) directly 
in the pair data type. 

d) When a search is performed on the pairs in A, B , ot Bj, full advantage 
can be taken of the fact that we may rely on data structures which allow 
logarithmic search costs. 

Looking at the timings above, we see that, on average and with comparable 
implementations, our new algorithm is faster than the Gebauer-Moller instal- 
lation. In some examples, the gains are relatively small, and in exceptional 
cases, the structure of the combinatorial data produces a larger overhead for 
our algorithm than for the Gebauer-Moller installation. 
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